Bragg spectroscopy of an accelerating condensate with solitary-wave behaviour 
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We present a theoretical treatment of Bragg spectroscopy of an accelerating condensate in a 
solitary-wave state. Our treatment is based on the Gross-Pitaevskii equation with an optical poten- 
tial representing the Bragg pulse and an additional external time-dependent potential generating 
the solitary-wave behaviour. By transforming to a frame translating with the condensate, we derive 
an approximate set of equations that can be readily solved to generate approximate Bragg spec- 
tra. Our analytic method is accurate within a well defined parameter regime and provides physical 
insight into the structure of the spectra. We illustrate our formalism using the example of Bragg 
spectroscopy of a condensate in a time-averaged orbiting potential trap. 



1. INTRODUCTION 

Bragg spectroscopy has proved to be a precise and ver- 
satile technique for measuring a variety of properties of a 
Bose-Einstein condensate, including the momentum dis- 
tribution [lj and coherence functions 2j. A theoretical 
proposal has been made to use Bragg spectroscopy to 
provide a signature of a vortex state in a condensate Q , 
and recently the technique has been suggested as a means 
for characterizing superfluid pairing in fermionic systems 
1. 

The theoretical description of Bragg spectroscopy of a 
condensate is well established and Blakie et al [j| , for ex- 
ample, have shown that a mean-field treatment using the 
Gross-Pitaevskii equation is an appropriate formalism. 
Unfortunately numerically calculating a Bragg spectrum 
using the Gross-Pitaevskii equation is highly computa- 
tionally intensive and may not readily yield a physical 
interpretation of the system. 

In this paper we present a treatment of Bragg spec- 
troscopy of a condensate which is accelerating in the form 
of a solitary- wave due to the influence of an external time- 
dependent trapping potential. Solitary-wave behaviour 
occurs only for a certain class of potentials, but is not 
uncommon. A simple example is the dipole oscillation of 
a condensate in a harmonic trap. Our method is based 
on transforming to the frame of reference where the con- 
densate is stationary, allowing us to establish a relatively 
simple set of equations that can be readily solved to give 
approximate Bragg spectra. We also obtain an analytic 
approximation to the spectrum, which provides physical 
insight into its structure. 

The paper is organised as follows. In Scc.[3]we describe 
our theoretical formalism, introducing solitary-wave so- 
lutions of the Gross-Pitaevskii equation and transforming 
the equation to the translating frame. We outline a sim- 
ple two-state momentum-space formalism which can be 
used to generate approximate Bragg spectra for a gen- 
eral solitary- wave, and we determine the validity regime 
in which such spectra are accurate. In Sec. [31 we use the 



methods of Sec. 0to treat Bragg spectroscopy of a con- 
densate undergoing micromotion in a time-averaged or- 
biting potential trap. We demonstrate the utility of our 
approximate solutions by comparing them with Bragg 
spectra calculated by simulating the Gross-Pitaevskii 
equation. Bragg spectroscopy in a time-averaged orbit- 
ing potential trap has been investigated experimentally 
|3 and our treatment provides quantitative agreement 
with full theoretical calculations and qualitative agree- 
ment with the experimental results 0. 



2. BRAGG SPECTROSCOPY OF A 
SOLITARY- WAVE 

Our theoretical treatment is based on solving the 
Gross-Pitaevskii equation for a condensate with wave 
function ip(r,t), in the presence of a time-dependent ex- 
ternal trapping potential and a Bragg pulse, i.e., 

■ r V>(r,t) - (-V 2 + V(r,t)+V opt (r,t)+ 



dt 



CVM)| 2 HM). 



(i) 



We have chosen to present our discussion in terms of di- 
mensionless units where the time scale is given by the 
inverse of a characteristic trap frequency lo x (e.g., as 
defined in Sec. 12. LIU and the position scale is xq — 
y^h/2muj x (where m is the mass of an atom) . The exter- 
nal trapping potential V(r, t) can be decomposed without 
loss of generality as 

V(r,t) = V(r)+G(r,t), (2) 

and the optical potential generated by the Bragg pulse is 



V op t(T,t) = -t/ cos(q-r- ujt), 



(3) 
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where q and u> are, respectively, the wave- vector differ- 
ence and the frequency difference between the two laser 
beams forming the Bragg pulse @- The dimensionless 
nonlinearity coefficient C is defined in terms of the num- 
ber of condensed atoms N, and the s-wave scattering 
length a, i.e., 

C = ^L. (4) 
mu) x XQ 
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2.1. Solitary- wave solutions 



2.2. Translating frame 



Solitary-wave solutions to the Gross-Pitaevskii equa- 
tion in the absence of the optical potential [Eq. with 
Uq = 0] exist for particular forms of the external po- 
tential V(r,t) of Eq. J5J. Such solutions, which evolve 
in three dimensions without changing shape, have been 
discussed previously 0, 0, |n| and have the form 



V-sw(r,i)=e(r-f(i)) e - vt+lS(r '* ) , 



(5) 



where the envelope wave function £(r) is an eigenstate of 
the time-independent Gross-Pitaevskii equation with the 
time- independent trapping potential V(r) and chemical 
potential /i, i.e., 

M C(r) = (-V 2 + y(r) + C|e(r)| 2 )e(r). (6) 
The position offset in the envelope wave function is 



^sw( r !*) r, 0sw (r,t)di 



J C(rH(r)dr. (7) 



The phase S(r, t) is determined by substituting the 
solitary- wave solution @ into the time-dependent Gross- 
Pitaevskii equation [TJ. We restrict our discussion to 
solutions where the velocity of the envelope £(r — f (t)) 
is proportional to the gradient of S(r, t) [in our units 
2V5(r, t)], and equating this to the centre of mass veloc- 
ity df(t)/dt yields 



S(r,i) 



dr(t) 
dt 



'(*) 



dt. (8) 



The centre of mass motion of the solitary-wave solutions 
can be shown to obey a form of Ehrenfest's theorem, 

1 d 2 f{t) 

2 dt 2 

where 



-VF(r,t), 



(9) 



F(r,t) = V(r,t) - V(r- f(t)). 



(10) 



For solitary-wave solutions to exist the function F(r, t) 
must yield a solution to Eq. © and, therefore, can be at 
most linear in r. 



2.1.1. Harmonic trap 

In Bose condensed systems it is most common that the 
time-independent trapping potential V(r) is harmonic, 



%(r) = i(* 2 + y 2 + AV), 



(11) 



with frequency lu x in the x-y plane (in SI units), aspect 
ratio A, and where r = {x, y, z). In the absence of the 
Bragg pulse (Uq = 0), the Gross-Pitaevskii equation JTJ 
has solitary-wave solutions provided that G(r, t) has the 
particular form such that Eq. © is satisfied, i.e., 



G(r,t)= Sl (t)-r + g 2 (t). 



(12) 



We begin our treatment of Bragg spectroscopy in a 
time-dependent external trapping potential by trans- 
forming to the frame translating such that the condensate 
centre of mass is stationary in that frame. While nu- 
merical solutions of the Gross-Pitaevskii equation can be 
calculated in the lab frame, the description of the system 
in the translating frame, where the condensate is at rest, 
allows simple approximate methods and a physical inter- 
pretation to be developed since the only time-dependence 
for the system in the translating frame appears via an ef- 
fective optical potential. 

A general quantum mechanical transformation to an 
accelerating translating frame has been given earlier |l l| . 
For the particular case we are interested in, of solitary- 
wave motion, the translation in co-ordinate space is de- 
fined by setting the coordinates of the new frame to be 



R = r-f(t), 



(13) 



where we denote the components of R to be (X,Y,Z). 
The momentum in the translating frame is derived by 
differentiating Eq. (|13fl yielding 



P = p 



1 df(t) 

2 dt ' 



(14) 



where we have used the fact that in our dimensionless 
units p = v/2. The wave function in the translating 
frame is 



^(R,*) = e-^^iPir^), 



(15) 



facilitating the momentum translation of Eq. i|14|) 17]. 
The Gross-Pitaevskii equation in the translating frame 
can now be derived by transforming Eq. using 
Eqs. (H2J and (O to yield 



ijypM = (-v 2 t + iy t (R,t) + v; t pt (R,t)- 



ci^CR,*)! 2 )^^,*)- 

In Eq. I|16|) we have denoted 



V 



R 



d d d 
dX'dY'dZ 



(16) 



(17) 



and the optical potential in the translating frame is 

l£ pt (R, t) = iu cos(q R + q f(t) - ut). (18) 

The potential denoted W t (R, t) is the sum of the external 
trapping potential in the translating frame and additional 
terms due to the momentum translation, i.e., 



VF l (R, t) = y(R + r(t),t) + -R 



1 d 2 r(t) 1 (dr(t) 



,ir 2 4 1" di 



(19) 
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By taking the gradient of Eq. I|19(l and making a com- 
parison with Eq. lO we find that 



V R W\R,t) = V R V(R) 



and, therefore, 



W^R, t) = V(H)+A(t), 



(20) 



(21) 



where A(t) is a function only of time. Defining an inter- 
action picture by 

ip T {K,t) = ^(R^tyS A{t)d \ (22) 
the Gross-Pitaevskii equation becomes 

;J^ T (R,<) = (-v^ + y(R) + y t pt (R,i)+ 

C|^ T (R,i)| 2 )^ T (R,i), (23) 

which shows that in the translating frame all the explicit 
time-dependence of the system is contained within the 
optical potential term [see Eq. (fT%f ]. When the Bragg 
pulse is turned off, stationary solutions of Eq. (|23|l exist 
of the form 



^ w (R,t) = £(R)e~ vt , 



(24) 



which arc solitary-wave solutions in the lab frame with 
centre of mass motion governed by Eq. @. 



2.3. Simple momentum formalism 

A theoretical formalism for describing Bragg spec- 
troscopy of a stationary condensate has been developed 
by Blakie et al. 0. Those authors derived an analytic 
treatment of Bragg scattering, for the case where C = 0, 
using an approximate two-state momentum space model 
to consider condensates released from a harmonic trap 
at the instant that the Bragg pulse is applied. Their 
treatment provides a simple approximate method for gen- 
erating Bragg spectra which, within a given parameter 
regime, are in quantitative agreement with full calcula- 
tions of the Gross-Pitaevskii equation. Here we extend 
that model to describe Bragg spectroscopy of condensates 
moving with solitary- wave behaviour. 

Following broadly the approach of Blakie et al. || we 
place three restrictions on the system. First, mean- field 
effects are neglected by taking C = 0. Our treatment of 
the nonlinear term is detailed in Sec. 12.41 Secondly, the 
time-independent trapping potential V(r) is neglected 
while retaining the time-dependent G(r, t). This means 
that the condensate centre of mass motion is still gov- 
erned by Eq. 10 but the condensate is no longer con- 
fined about its centre of mass. When neglecting atomic 
interactions, removing the trap V(r) does not alter the 
momentum distribution of the unscattered condensate, 
whereas the condensate part scattered by the Bragg pulse 
is decelerated. Provided that the Bragg pulse length is 



short compared to the oscillation period of atoms con- 
fined by V(r), the effect of the confining potential on 
the Bragg spectrum is negligible. In some experimen- 
tal systems it may be possible to apply V(r) and G(r, t) 
independently so that V(r) can be turned off when the 
Bragg pulse is applied. Thirdly, for simplicity, we choose 
the wave vector q of the optical potential to lie along the 
x axis such that q • R = qX and the optical potential of 
Eq. i|18|) becomes 



-U cos(qX + qx(t) - cot). (25) 



This can be done without loss of generality for the case 
where V(r) is turned off. Finally then, within these re- 
strictions, the wave function i[i T (R, t) is seperable and 
denoting the product wave function in X by ip T (X, t) it 
is possible to reduce the analysis of the system to one 
dimension without loss of generality Q • 

We begin by transforming the one-dimensional equiv- 
alent of Eq. J2S1), with C = and V(R) = 0, into mo- 
mentum space using the Fourier integral 



cf> T (K,t) 



'2it 



^ T (X,t)e- 



iKX 



dX, 



(26) 



to give 

a 



e-i^t-im)^^ - q,t) 



+e %{ujt - qS{t)) 4> T {K + q,t) 



(27) 



Next, we partition momentum space into bins of width 
q, centred around the Bragg orders K = nq, where n is 
integer. That is achieved by expressing the wave function 
4> T (K, t) as a set of wave functions 4>n(K, t), each defined 
only within bin n of the partitioned momentum space, 



4%(K,t) =f(K,t), (28) 



where 



K;(n-^jq<K<(n+^\q. (29) 

The description of the system is simplified by defining 
k = K — nq and writing 



in(ujt — qx(t)) 



(30) 



where k varies from —q/2 to q/2, for every bin. The mo- 
mentum space Schrodinger equation can then be written 
as an infinite set of equations coupling wave functions of 
consecutive bins, 

d 

+^M$-iM) + 0JE+i(*.*)] >( 31 ) 
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where 



u) n (K, t) — (k + nq) 2 — nuj + nq 



dx(t) 
dt 



(32) 



Equation is identical to that derived by Blakie et 
al. for a stationary untrapped condensate, with the 
exception of an additional time-dependence appearing 
in oj n (K,t). The additional term nqdx(t)/dt is due to 
the momentum of the condensate in the lab frame [see 
Eq. and has previously been referred to as the 

Doppler term [6j. 

To determine the Bragg spectrum of an accelerating 
condensate we calculate the population in the first Bragg 
order (n = +1 bin) as a function of the frequency differ- 
ence lo in the optical potential. In the translating frame 
this is given by 



P+i(w,t) 



corresponding in the lab frame to 

P+l(w,t)= / \cb(k,t)\ 2 dk, 



where <^(fc,i) is defined by 
1 



ip(x,t) e - ikx dx 



(33) 



(34) 



(35) 



in terms of the wave function in the lab frame tp(x,t), 
and the limits of integration in Eq. Ij34(l account for the 
solitary-wave centre of mass motion of the condensate. 

Solving Eq. l(3*T|) . in order to calculate (/^(k, i), is no 
more straight-forward than solving the full Schrodinger 
equation. Again following the treatment of Blakie et 
al. Q the problem becomes readily tractable if we assume 
that in the region of interest, where first order Bragg 
scattering is strong, only the n = and n = +1 momen- 
tum wave functions are substantially populated. In this 
'two-state model' Eq. (|31|l reduces to 



d 



0o (X t) 

4>+i (M) 



rU 



jU u> (k) + A ^+i(k, t) 



(f>o(K,t) 

0li(M) 



(36) 



where loq(k) = k 2 and Ao^+i(k, t) = q 2 + 2nq — lu + 
qdx(t)/dt. Equation (|36|l can be solved easily by numer- 
ical means and its solution accurately describes the sys- 
tem when n ^ 0, +1 populations are negligible. The 
validity conditions for Eq. I|36(l are now derived. 



2.3.1. Validity of the two-state model 



for two-photon transitions only (which couple consecu- 
tive bins) the detuning for the transition from bin n to 
bin m is defined by 



A, 



t (K, t) = LO m (K, t) - CJ„(K, t), 



(37) 



where \m — n\ = 1. We define a transition to be on 
resonance when the magnitude of the detuning is less 
than Uq/2 (determined by the power-broadened width 
@). Thus, the first order Bragg transition (between the 
n = and n = +1 bins) will be resonant for some part 
of the initial condensate momentum distribution over the 
frequency range 



and 



, dx(t) 

UJ>q+q ^r 



n dx(t) 

lu <q 2 + q 



1 



dt 



aq 



aq 



:U 



(38) 



(39) 



where a is the full momentum width of the condensate. 
We have chosen to use the 1/e condensate momentum 
width and the extreme values of dx(t) / dt in order to yield 
the largest possible transition width and, therefore, the 
strictest condition for the two-state validity. Within the 
frequency range defined by Eqs. (|38|l and (|39|l we require 
that the transitions to neighbouring bins (n = — 1 and 
n = +2) are not resonant, i.e., |A_i„» ( K , 01 > Uq/2 and 
|A+i_>+2(/t, t)\ > C/o/2. Those conditions are satisfied 
when 



U < 2q(q -a - 5 P ), 



(40) 



where 25 p = dx(t)/dt\ max — dx(t) / dt\ m - ln . Note that 
Eq. i|4U|) requires a < q, which also insures that momen- 
tum wave functions of consecutive bins do not overlap. 

The discussion leading to Eq. (|40|l does not account 
for higher order processes. It is possible for 2m— photon 
processes to transfer population from the n = to the 
n = m bin, without intermediate bins accumulating sig- 
nificant population. However, provided the Rabi period 
for such processes is significantly longer than the Bragg 
pulse length, the population transfer is negligible and will 
not effect the two-state dynamics between the n = and 
n = +1 bins. From the higher order processes, the four- 
photon processes are the most likely to effect the two- 
state validity. In the case of no centre of mass motion 
(x(t) — for all t), the Rabi frequency at the resonance 
of the n = to n = +2 transition is 



Rabi(0^+2) 



If 2 



II 2 



8(u/ (0)-u>+i(0)) 



(41) 



Complete population transfer (for n = 0) would occur for 
^Rabi(o^+2)^ = 71 'j yielding the additional condition 



The two-state model is valid when only the n = and 
n = +1 bins are significantly populated. Accounting 



(42) 



5 



The condition derived equivalently on the n = — 1 to 
n = +1 transition is more lenient and therefore, provided 
that Eqs. JUJ and JUJ are satisfied, all other transitions 
can be neglected and the two-state approach is valid. 



2.3.2. Simplified treatment of the two-state model: Main 
features of the Bragg spectrum 



The strength, relative phase, and resonance position for 
each term in the sum can be determined by finding q, 
ei, and u>i, respectively. Those parameters can be calcu- 
lated relatively easily by equating Eqs. 1251) and Q48[). de- 
composing the former into trigonometric functions, and 
comparing terms. This leads to 



iuiit+iei 



(49) 



The two-state model provides approximate Bragg spec- 
tra for a condensate with solitary-wave behaviour and 
requires much less computational resources than calcu- 
lations of the full Gross-Pitaevskii equation. We will 
demonstrate the success of the two-state solutions to 
Eq. i|36[l m Sec. 13.31 However, the dominant features 
of the Bragg spectrum can be readily obtained from a 
particular simplification of the two-state model, which 
we outline in this section. 

Blakie et al. ||| have shown that if the optical poten- 
tial (|18l) has the simple form 6cos(q • R — u)t)/2, then 
for the case F(R) = and C = 0, the detuning in the 
two-state model becomes time- independent, i.e., 



A ^ +1 (k) = q 2 + 2nq - u. 



(43) 



In that case, the analytic solution to Eq. (|36(l yields the 
Bragg spectrum 



y 6 2 |0 O ( K ,O)| 2 



where the Rabi frequency is 



[1 - COS f^Rabi («)*] dn, 



^Rabi(«) 




+ A2 



(44) 



(45) 



The frequency width of the first order Bragg transition 
is 



= (- 



(46) 



where contributions from power broadening and tempo- 
ral broadening are included. In the limit bt/4 <C 1 the 
total scattered population across the single Bragg reso- 
nance is 



P+i (lj, t)duj 



nbH 



(47) 



For the more general case we are considering here, the 
optical potential in the translating frame is given in one 
dimension by Eq. (|25|l . but it can be expanded as a linear 
sum of simple Bragg potentials with the form considered 
by Blakie et al. @, i.e., 



^^2ciaoa(qX-(u + ui)t + ei). (48) 



and taking the fourier transform of both sides of Eq. I|49|) 
we find that 



U, 



/CO 
e iqx(t) e -iwt dt = J2 cie ltl 8{w + ui), (50) 
-OO j 



which allows q, e/, and lui to be readily extracted, for 
all I. 

The Bragg spectrum of a condensate accelerating as 
a solitary-wave will have a number of resonance peaks 
due to the multiple Bragg terms in the linear sum of 
Eq. H48fl . When the individual Bragg spectra due to 
consecutive terms in frequency space overlap, the effect 
due to simultaneous application of the Bragg fields can 
be complicated. However, there are particular regimes 
where the Bragg response at each ui can be considered 
as if it were acting independently of the others, and then 
Eqs. (|43|) to (|47|l can be applied to each resonance, and 
it is straight-forward to construct a complete approxi- 
mate spectrum by summing the Bragg spectra due to 
each individual term. We refer to this method of calcu- 
lating approximate spectra as the independent resonance 
approach. The most important case where this approach 
applies is when the position of consecutive resonances are 
separated by more than the widths of those resonances, 
i.e., 



LU l + 1 ~U>1> 



er w (c;+i) + cr^q) 



(51) 



Other particular cases arise, for example in the perturba- 
tive limit q£/4 <§; 1, the independent resonance approach 
is accurate when 

cos -(oJi+i - u}[)t - (ej+i - ej) =0. (52) 



2.4. Nonlinear effects 

Calculating numerical solutions to the three- 
dimensional time-dependent Gross-Pitaevskii equation is 
computationally intensive. Blakie et al. @ demonstrated 
that within a well defined parameter regime the analytic 
results derived for the case where C — provide a good 
representation of the main behaviour observed in the 
many-body system. For the systems we are considering, 
where the condensate has solitary-wave behaviour in the 
absence of the Bragg pulse, the centre of mass motion of 
the condensate is independent of the nonlinear strength 



6 



of the system and, therefore, in the translating frame 
the validity regime definitions derived by Blakie et al. Q 
apply. 

There are three dominant many-body effects on Bragg 
scattering. The first is that atomic interactions within 
the condensate cause the momentum distribution to 
broaden during free expansion. The time scale for this 
has been derived to be t = 7%/3/4/x [8j and provided that 
the Bragg pulse length is shorter than this time scale, 
the effect of expansion is negligible. The second many- 
body effect is a shift in the resonance condition of the 
first Bragg order due to the mean-field. For condensates 
confined in a three-dimensional harmonic trap the reso- 
nant frequency is higher by approximately 4/x/7 @ and 
this shift can be neglected when /i -C Uq/2. In addition 
to these effects, harmonically trapped condensates with 
a higher chemical potential have a narrower momentum 
distribution. This arises due to the form of £(R), de- 
fined by Eq. ©: a broader spatial distribution leads to 
a narrower momentum distribution. The Bragg spectra 
for trapped condensates with a higher chemical poten- 
tial exhibit narrower spectral features and increased peak 
scattering. These effects were also observed by Blakie et 
al. 0. 

In the case where C — a description of the three- 
dimensional system can often be simplified to one dimen- 
sion without loss of generality, as described in Sec. 12.31 
The Gross-Pitaevskii equation, while it is not separa- 
ble, can also be reduced to a one-dimensional form (e.g., 
[l3L Pill ID^I. Bragg spectra calculated using a one- 
dimensional treatment along the axis of the Bragg pulse 
also generates reasonably accurate results for a nonlin- 
ear condensate, as we will see for a particular example in 

Sec. cm 



3. TIME-AVERAGED ORBITING POTENTIAL 
TRAP 

In this section we calculate Bragg spectra of a con- 
densate accelerating as a solitary-wave. The spectra are 
calculated at various levels of approximation: (i) using 
the approximate methods described in section Sec. 12.31 
(ii) using the full Schrodinger equation (Gross-Pitaevskii 
equation with C = 0), and (iii) using the Gross-Pitaevskii 
equation in both one and two dimensions. 

Condensates confined by a harmonic trap will exhibit 
solitary-wave behaviour for an external potential satisfy- 
ing Eq. (|12p. Some examples include the dipole oscilla- 
tion of a condensate in a harmonic trap, and the motion 
resulting from translations of the trap due to the effect of 
noise. We have chosen the solitary-wave micromotion of a 
condensate in a time-averaged orbiting potential (TOP) 
trap to illustrate our methods. In this particular system 
the condensate moves in a circular trajectory following 
the bias field that removes the discontinuity in the trap- 
ping potential (Jill and references therein). The micro- 
motion velocity of the solitary-wave is so high that the 



structure of the Bragg spectrum can be measured exper- 
imentally with moderate frequency resolution. The ex- 
perimental investigation of Bragg spectra in a TOP trap 
|6| shows qualitative agreement between the experimen- 
tal data and calculations made using the one-dimensional 
Schrodinger equation. The application of our methods of 
Scc.|21provide the theoretical framework for the concepts 
introduced previously in relation to Bragg spectroscopy 
in a TOP trap HHg. 

3.1. Theoretical description of the TOP trap 

The TOP trap potential can be well approximated by 

V^top (r, t) = V H (r) + r (cos fit, sin fit, 0) • r, (53) 

where r$ is the so-called 'circle of death' and fi is the 
bias field rotation frequency Equation (|53|l is of 

the form required for solitary-wave solutions to exist [see 
Eq. (|12[l] and in fact condensation occurs into the sys- 
tem ground state, which is a solitary-wave solution of 
the Gross-Pitaevskii equation with circular centre of mass 
motion given by 

f(t) = 7t (cos fit, sin fit, 0), (54) 

where 7t = 2r /(fi 2 - 1) 

3.2. Effective optical potential 

The solitary-wave description of condensate micromo- 
tion allows Bragg spectroscopy in a TOP trap to be de- 
scribed using the theory presented in Sec. [3 In the trans- 
lating frame the optical potential of the Bragg pulse be- 
comes, in one dimension, 

V* pt (X,t) = -U cos(qX + q"f t cosnt-u;t). (55) 

The optical potential of the Bragg pulse, in the moving 
frame of the condensate, has been identified prev iously 
as being a frequency modulated potential 0, 161. Our 
discussion in the translating frame provides a theoretical 
validation of that identification. 

In anticipation of interpreting Bragg spectra of a con- 
densate in a TOP trap we decompose the optical poten- 
tial in the translating frame in the form of Eq. (|48Jl . For 
the frequency modulated potential of Eq. (|55|l such an 
expansion is well known, i.e., 

V* pt (X,t) = ~U {J { lq )co S (qX-wt) 

OO 

- ^(-1)' J2/+1 (79) [cos (qX - vt 
1=0 

-(2Z + l)fit- 

+ cos (qX - Lot + (21 + l)fif - D 
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FIG. 1: Bragg spectrum of a condensate in a TOP trap 
calculated using the full Schrodinger equation (-), the two- 
state model (- -), and the independent resonance approach 
(■■). The vertical lines along the top of the figure indicate the 
positions of the spectral resonances predicted by Eq. 1561 , 
Parameters are Uo — 45, q — 29, ro = 1241, £1 = 153, and 
t = 0.105tt (8 TOP trap cycles). 




FIG. 2: Bragg spectrum of a condensate in a TOP trap 
calculated using the full Schrodinger equation (-), the two- 
state model (- -), and the independent resonance approach 
(••). The vertical lines along the top of the figure indicate 
the positions of the spectral resonances predited by Eq. 1561 1. 
Parameters are Uo — 45, q = 29, ro = 878, Q = 76, and 
t = 0.092tt (3.5 TOP trap cycles). 



+ ^(-1)' J 2 ;(7g) [cos(qX -uxt- 2lttt) 
1=1 

+ cos{qX - ujt + 2mt)}} , (56) 

where J u (z) is the Bessel function of the first kind of or- 
der v and with argument z. In the translating frame, 
the condensate is subjected to an optical field which is 
the sum of simple Bragg terms, separated in frequency by 
the TOP trap bias field rotation frequency. Therefore, we 
expect that the Bragg spectra of a condensate with mi- 
cromotion will have multiple resonances separated by fl. 



3.3. Bragg spectra 

The Bragg spectrum of a condensate in a TOP trap is 
presented in Fig. ^ for a parameter set which is exper- 
imentally accessible p. The spectrum calculated using 
the full Schrodinger equation (i.e., C = 0) is shown as 
a solid line. The Bragg spectrum is centred at the first 
order Bragg resonance [see Eqs. (|38|) and (|39[l ] and ex- 
hibits additional structure compared to the single reso- 
nance peak at u> = q 2 which would be expected from a 
single Bragg term. Equation (|56[) correctly predicts the 
positions of the resonance peaks in the Bragg spectrum 
(as indicated by the vertical lines in the figure). 

We have also generated approximate Bragg spectra of 
a condensate in a TOP trap using the two-state model 
(see Sec. 12. 3|) . as well as the simplified version of the 
two-state model, the independent resonance approach 
(see Sec. 12. 3.211 . Solving the two-state equation i|3T)|) 
is two orders of magnitude faster than solving the full 



Schrodinger equation. Provided that the validity condi- 
tions of Eqs. H4U|) and (|42|l are satisfied, the results of 
the two-state model are reasonably accurate, as demon- 
strated by the dashed line in Fig. ^ The main discrep- 
ancy between the two-state method and the full C = 
calculation is the underestimation of the two strongest 
resonance peaks. This occurs because the stationary har- 
monic trap in the translating frame has been neglected 
in the two-state approach. The effect of the harmonic 
trap is to shift population towards zero momentum thus 
allowing more condensate population to be scattered into 
the n — +1 momentum order. The calculation using the 
independent resonance approach is indicated by the dot- 
ted line in Fig. ^ and is almost identical to the Bragg 
spectra calculated using the two-state approach. This is 
because each Bragg resonance within the spectrum is well 
resolved [refer to Eq. H51[)] and, futhermore, the spectrum 
in Fig.^is for exactly eight TOP trap cycles and we note 
that Eq. (|52|l holds in that case for TOP trap resonances 
which are consecutive in frequency space. 

Figure [5] shows a Bragg spectrum for a condensate in 
a TOP trap where the bias field rotation frequency is 
halved compared with the case of Fig. ^ and the Bragg 
pulse length is not an integer multiple of that frequency. 
The spectrum is again centred at uj — q and the reso- 
nances are separated by the TOP trap rotation frequency, 
although a number of peaks appear missing (including 
the central one) because their weight is small in expan- 
sion (|56|l . In contrast to Fig. ^ the Bragg spectrum of 
Fig- EJis no longer symmetric about the centre of the first 
order Bragg resonance, but favours the low frequency res- 
onances due to our initial conditions which yield a neg- 
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FIG. 3: Bragg spectrum of a condensate in a TOP trap 
with the parameters of Fig.^calculated using the Schrodinger 
equation (-), the one-dimensional Gross-Pitaevskii equation 
with Cm = 85, n = 10.08 (- -), and the two-dimensional 
Gross-Pitaevskii equation with C2D = 600, fi — 9.85 (■•). The 
vertical lines along the top of the figure indicate the positions 
of the spectral resonances predicted by Eq. H56I . 



ative Doppler term for the first half of the TOP trap 
rotation. We observe that the two-state model does ex- 
hibit the lack of symmetry present in the full Schrodinger 
equation calculation, while the independent Bragg reso- 
nance approach does not. We also note that the two-state 
model more accurately approximates the full Schrodinger 
equation calculation than was the case in Fig-QJ The rea- 
son for this is that neglecting the harmonic trap has less 
effect at lower scattered fractions. 



3.4. Effect of condensate nonlinearity 

In this section we present Bragg spectra for a nonlinear 
condensate in a TOP trap calculated using the Gross- 
Pitaevskii equation in both one and two dimensions. In 
general the Gross-Pitaevskii equation is not separable so 
we determine the appropriate nonlinear strength in the 
reduced dimensional system by arranging to have the 
same Thomas-Fermi radius and chemical potential as in 
the three-dimensional description. 

Figure [3] illustrates the Bragg spectra for the same 
parameters as in Fig. ^ but with C ^ 0. The nonlin- 
ear strength corresponds to approximately 2 x 10 4 atoms 
in the Otago TOP trap 0. The results from the one- 
dimensional Gross-Pitaevskii equation agree well with 
the Bragg spectrum calculated in two-dimensions, jus- 



tifying our simplification of the problem to one dimen- 
sion. Calculating the two-dimensional Bragg spectra in 
Fig. n is computationally expensive, requiring approxi- 
mately 10 CPU hours on a 3 GHz Pentium to calculate 
each point (with 100 points required for the full spec- 
trum). 

We have also reproduced in Fig. U| the result from the 
linear case of Fig. ^ an d comparing the C — and the 
C / spectra the effects discussed in Sec. 12.41 are ob- 
served. The nonlinear Bragg spectra exhibit increased 
peak scattered fractions, reduced resonance width, and a 
slight increase in the resonant frequencies, with respect 
to the equivalent C = system. The increased resolu- 
tion in the features of the Bragg spectrum occurs due 
to the reduced momentum width of the initial state, de- 
termined by the form of £(R). It is interesting to note 
that the two-state model can also exhibit the narrowed 
resonances if the initial state is taken to be the same nu- 
merically calculated ground state as was used for the full 
nonlinear Gross-Pitaevskii calculation. 



4. CONCLUSION 

We have presented methods for identifying the key fea- 
tures of the Bragg spectrum of a condensate which ex- 
hibits solitary-wave behaviour. In particular such an ac- 
celerating condensate becomes stationary in a translating 
frame and an effective optical potential contains the only 
explicit time-dependence of the system in that frame. We 
have generalized the two-state formalism introduced by 
Blakie et al. @ and have defined the regime in which the 
approximate spectra generated using that approach are 
in good agreement with full calculations. Bragg spec- 
tra of a solitary- wave exhibit multiple resonances and we 
have identified particular cases where those resonances 
each act independently of the others, thereby revealing 
the main features of the Bragg spectrum with minimal 
calculation. 

The Bragg spectrum of a condensate in a TOP trap 
has been used to illustrate the application of our the- 
ory. We have described the resonance structure in the 
spectra and in the validity regime of our treatment have 
demonstrated excellent agreement between our approxi- 
mate methods and full theoretical calculations. 
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